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ABSTRACT 

Aims. In this paper we study the effects of a toroidal magnetic flux tube emerging into a magnetized corona, with an emphasis on 
large-scale eruptions. The orientation of the fields is such that the two flux systems are almost antiparallel when they meet. 
Methods. We follow the dynamic evolution of the system by solving the 3D MHD equations using a Lagrangian remap scheme. 
Results. Multiple eruptions are found to occur. The physics of the trigger mechanisms are discussed and related to well-known 
eruption models. 
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1. Introduction 

The trigger mechanism for large-scale eruptions, such as coro- 
nal mass ejections (CMEs), is one of the main concerns of 
theoretical solar physics. There are many different approaches 
to modelling solar eruptions, ranging from 2D analytical mod- 
els to 3D n umerical simul ations. A review of such models can 
be found in lForbes| (1200 0) and further references are given in 
MacTaggart & Hood (2009a). One high ly influential model fo r 
solar eruptions is the breakout model (lAntiochos et al.l fl999l) . 
This theory uses, as an initial condition, a multipolar configura- 
tion, which consists of four distinct flux systems separated by a 
null point. This configuration is stressed by an imposed shearing 
at the lower boundary (normally taken to be the photosphere). 
Reconnection across the current sheet at a stressed null changes 
the field geometry and weakens the tension of overlying field 
lines, converting them into nonrestraining field lines of neigh- 
bouring flux systems. The continuation of this process can even- 
tually lead to expulsion of a flux rope. For more details of th e 
breakout model, the reader is directed to Antiochos et aL (1999), 
iDeVore & Antiochosl (l2008h andlLyiich et al.l (120081). 

Recently, IDeVore & Antiochosl (120081) and iLynch et all 
(2008) have simulated magnetic breakout in 3D. The first of 
these studies finds homologous confined eruptions. The second 
study follows the to pological evolution of a fast breakout CME. 
ISoenen et al.l ([2009) simulate homologous CMEs in the solar 
wind in an axisymmetric 2.5D configuration. 

The breakout studies mentioned above have been very useful 
in investigating the physical mechanisms of solar eruptions. One 
point which they all have in common, however, is that the initial 
equilibrium is stressed by artificially imposed shearing motions. 

Other models are also used to study eruptions in the solar 
atmosphere. For example, Mackay & van Ballegooijen (2006) 
model the large-scale coronal field as a series of non-linear force 
force-free equilibria. Flux ropes form but not all of them set- 
tle into equilibrium. Those that diverge from equilibrium are 
ejected. 

For studying dynamic evolution in the atmosphere, another 
model is flux emergence. This considers the early evolution of 
active regions, which, of course, are the sources of large-scale 
solar eruptions. The standard practice in dynamic flux emer- 



gence experiments is to have a stratified solar atmosphere in- 
cluding the top of the solar interior. A flux tube is placed in 
the solar interior and the system is left to evolve by itself, with 
no imposed flow s. A review of such models can be found in 
lArchontisl (f2008). Some studies inc lude a magnetized co r ona in 
their model. Archon tis et al.l (120051) and|Galsgaard et al. (2007) 
study the effects of a cylindrical flux tube emerging into a hori- 
zontal coronal field. They consider different orientations for the 
magnetic fields and discuss the re connection and high-speed jets 
that occur. iMaclean et al.l (120091) consider the same experiment 
but investigate the magnetic topology of th e system. In relation 
to solar eruptions, [Manchester et al.l (|2004) report on the erup- 
tion of a flux rope that forms during the eme rgence of a cylindri- 
cal flu x tube into an non-magnetized corona. lArchontis & Torokj 
(120081) study the emergence of a cylindrical flux tube into a 
magnetized corona and find that with a favourable orientation, 
a CME-like eruption is possible. 

In the present paper, we shall consider a similar setup to 
lArchontis & Torok (2008) but use a different model for the mag- 
netic field. This is a toroidal loop, placed in the solar interior, 
rather than a cylindrical one dHoodet_am009l A comparison o f 
these two models is discussed in iMacTaggart & Hoodl (H)09b). 
This paper will show that multiple large-scale eruptions are pos- 
sible from the same emerging region and will discuss the physics 
of the trigger mechanisms. 

The outline of the paper is as follows: §2 will describe the 
model setup and the initial conditions. In §3 we shall discuss 
the results of the eruption experiment and link the processes in- 
volved with eruption models. §4 will summarize the results. 



2. Model setup 

To study the evolution of the system num erically, we use a 
Lagrangian remap scheme (lArber et al.ll200ll) . All variables are 
non-dimensionalized with photospheric values: pressure, p V h = 
1.4 x 10 4 Pa; density, p p h = 3 x 10~ 4 kg irT 3 ; temperature, 
r ph = 5.6 x 10 3 K and scale height H ph = 170 km. The other 
units used in the simulations are: time, f p h = 25 s; speed, u v h = 
(/?ph/p p h) 1/2 = 6.8 km i" 1 and magnetic field B p h = 1.3 xlO 3 G. 
The evolution of the system is governed by the following resis- 
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tive and compressive 3D magnetohydrodynamic (MHD) dimen- 
sionless equations 

dp 

p{^ +(u- V)uJ = -V/? + (VxB)xB+pg, 



m 

dt 



V x (u x B) + 77V 2 B, 



P |^ + (U ■ V)ej = -pV ■ u + T]f, 
V ■ B = 0, 

with specific energy density e = p/((y - l)p). The basic vari- 
ables are the density p, the pressure p, the magnetic field vector 
B and the velocity vector u. j is the magnitude of current den- 
sity and g is gravity (uniform in the z-direction). y is the ratio 
of specific heats and is taken as 5/3. r\ is the resistivity which is 
taken to be uniform with a value of 0.001. The code accurately 
resolves shocks by using a combination of artificial viscosity and 
Van Leer flux limiters. In such regions, heating is added to the 
energy equation. 

The equations are solved in a Cartesian computational box 
of (non-dimensional) sizes, [-50,50] x [-50,50] x [-20,85] 
in the x,y and z directions, respectively. The boundary condi- 
tions are closed on the top and base of the box and periodic 
on the sides. The background stratification is similar to that in 
IMacTaggart & Hoodl d2009b). The non-magnetic stratification 
includes a solar interior that is marginally stable to convection 
(-20 < z < 0), a photosphere/chromosphere (0 < z < 10), a 
transition region (10 < z < 20) and a corona (20 < z < 85). A 
horizontal magnetic field of the form 

B = B c (z)( 1,0,0), 

is included in the corona. B c {z) is a hyperbolic tangent profile, 
so that the field is uniform in the corona and rapidly declines 
to zero at the base of the transition region. The strength of the 
coronal field is taken to be 0.01 (« 13G). The orientation of the 
field is chosen so that it is almost antiparallel to the field of the 
emerging flux tube when they meet. The initial toroidal tube, that 
is placed in the solar interior, has the form 



B x = B g (r) 



s - s 



B v = -B^r) Z -^-B e (r)- y -, 
s r s 

B z = B^-Betf)*^, 
s r s 

with 



x + (s - so) , s - so — rcosfl, x = rsin6>, 



and 



Ba - B e 



B(f — arBj, — aB§re 



so is the major axis of the tube and ro is the minor axis, zo is 
the base of the computational box. a is the initial twist and Bo 
is the initial axial field strength. Varying these parameters can 
have a profound effect on th e behaviour of toroidal emergence 
dMacTaggart & Hoo d 2009b). For this study, however, we shall 
only consider the values of Bo = 5 and a = 0.4. In this study 
we take so = 15, ro - 2.5 and zo = 25. To initiate the experi- 
ments, the entire tube is made buoyant, i.e. a density deficit rel- 
ative to the background den sity is introduced dHood et al.ll2009t 
iMacTaggart & Hood 20091). 



3. Results 

As mentioned previously, the pre-existing field of the corona is 
almost antiparallel to the field of the emerging tube. As the tube 
emerges, the arcade it forms first makes contact with the coronal 
field at t « 33. It pushes into the horizontal field and an arched 
current sheet forms between them. Reconnection occurs and the 
outer field lines of the arcade, that pass through the current sheet, 
change their connectivity and connect to the coronal field. The 
form of reconnectio n is fully 3D in the s ense that it does not 
involve a null point. iMaclean et alj d2009) also find this to be 
the case, with two clusters of null points forming at the sides 
of the emerging arcade rather than at the apex. They suggest 
that separator reconnection occurs in the current sheet between 
the tube and the corona. The effect, however, of this external 
reconnection is to weaken the tension of the coronal field. This 
becomes increasingly rapi d with the loss of t he restraining field, 
as in the breakout model dLvnch et al.ll2008l) . Evidence for this 
can be found by looking at the reconnection rate as it varies in 
time. The rate of reconnected flux is given by 



dt 



-I 



E«dl, 



where the right hand side is the integrated pa rallel electric field 
along the reconnection line (ISchindlerj [2008). In 3D, however, 
no unique fine exists a t which the flux is split and reconnected 
dHornig & Priestll2003l) . Within the current sheet, every field line 
constantly changes its connection. Therefore, to estimate the re- 
connection rate between the flux tube arcade and the corona, in 
the simulation, we measure 



dt 



max(£ , |i)(5, 

v=0 



where max v= o(£'||) is the maximum parallel electric field in the 
y = plane (at the top of the arched current sheet) and 5 is the 
thickness of the current sheet. This is a conservative estimate of 
the reconnection rate, giving a lower bound. Figure [TJ displays 
the increasing external reconnection rate d uring the expansion 
of the emerging arcade. iLvnch et afl (|2008) report that this be- 
haviour is also found in the breakout model. Figure [2] shows the 
field line structure of the reconnected corona at t — 70. The 
coronal field lines are traced from opposite sides of the com- 
putational box (red for one side (x = -50), green for the other 
(x = 50)) at z — 22. A magnetogram showing B- is placed at the 
base of the photosphere (z = 0). After external reconnection the 
coronal field connects down into the main photospheric polari- 
ties (sunspots), leaving the centre free for the arcade to push on 
upwards. 

As the arcade emerges, the Lorentz force produces shear- 
ing along the polarity inversion line between the two sunspots. 
Also during the rise of the arcade, plasma drains down the field 
lines, with some flowing into a region of reduced pressure within 
the arcade. This, combined with the shearing produced by the 
Lorentz force, results in internal reconnection in the arcade and 
the production of a flux rope. In the simulation this occurs at 
t « 80. For more details of the physics of the internal re- 
connection and flux rope pr oduction, the reader is directed t o 
IMacTaggart & Hoodl d2009bh (see also lManchester et allllool . 

With the production of a flux rope and the continued loss 
of the restraining coronal field through external reconnection, 
there is a catastrophic expulsion of the flux rope at t « 91. This 
description of the eruption is very similar to that of magnetic 
breakout. There are two important differences, however. The first 
is that in this flux emergence model, a null point does not play a 
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Fig. 1. As the arcade of the emerging tube expands into the mag- 
netized corona, the continued weakening of the overlying re- 
straining field results in faster reconnection. 




Fig. 2. The field line structure of the corona after external recon- 
nection. Red and green coronal lines are traced from opposite 
sides of the box, x = -50 and x = 50, respectively, at z — 22 and 
t = 70. A magnetogram shows B z at z — 0. 



crucial role in the external reconnection, as it does in breakout. 
The second is that there are no imposed flows in the flux emer- 
gence model. The Lorentz force of the emerging arcade naturally 
shears the magnetic field (MacTagg art & Hoo d 2009b). 

After the first eruption, the system does not settle into an 
equilibrium since the flux tube continues to emerge. The mag- 
netic field in the corona is no longer horizontal due to recon- 
nection during the first eruption. Shearing in the emerging ar- 
cade produces a second flux rope at a height z ~ 34 and time 
t as 125. The field line structure of the new flux rope and the 
corona at t — 125 is displayed in Figure [3] Coronal field lines 
(cyan) are traced from opposite sides of the box at z = 40. The 
red field line is the axis of the new flux rope. Some surround- 
ing field lines (green) are traced from the main photospheric 
polarities, which are shown on a magnetogram at z = 0. This 
flux rope also erupts but the mechanism of the eruption is dif- 




Fig. 3. The field line configuration just before the eruption of the 
second rope. Coronal field lines (cyan) are traced from the sides 
of the box (along the y-axis) at z = 40. The red field line is the 
axis of the newly formed flux rope. Some field lines surrounding 
this are traced in green. A magnetogram is displayed at z = 0. 



ferent to that of the first. The first eruption was of breakout type 
where significant external reconnection (between the emerging 
arcade and the horizontal coronal field) played a crucial role 
in the expansion phase up to the eruption. As can be seen in 
Figure [3] however, the first eruption has 'cleared a path' in the 
corona. External reconnection, therefore, does not feature as a 
key factor in the eruption of the second rope. This means that 
the trigger mechanism is no t of breakout type and othe r candi- 
dates need to be considered. Archontis & Torok (2008) suggest 
that the eruption they found is mainly driven by runaway recon- 
nection below the flux rope. This process has two important ef- 
fects on the eruption. Firstly, the reconnection weakens the ten- 
sion of the overlying field. Secondly, the upward reconnection 
jet carries the reconnected field li nes to the eru pting flux rope, 
adding poloidal flux to the rope (Vrsnak 2008). Both of these 
effects help the rope to accelerate further. This, in turn, enhances 
the reconnection below and a runaway process ensues. For the 
second flux rope, this in combination with a weakened coronal 
field could drive the eruption. Another possibility, however, is 
that the newly formed rope becomes subject to an ideal MHD 
instability. A likely candida t e is the torus instabili ty dBatemanl 
Il978t iKliem & T6r8gi2006t iTorok & Klie"rn1r2007h . A toroidal 
flux tube will become unstable against expansion if the exter- 
nal poloidal field decreas es su f ficient ly rapidly in the direction 
of the major tube radius. iFanl d2009l) simulates the emergence 
of a cylindrical flux tube into a non-magnetized corona. The au- 
thor claims that the expansion and acceleration of the flux rope, 
produced in the simulation, is due to the continuous injection of 
twist, from the interior, via torsional Alfven waves. 

Identifying the trigger mechanism for the second eruption is 
not a trivial undertaking and is therefore beyond the scope of 
this research note. This is a subject for future study. It is possible 
that there is no single trigger mechanism and that the second 
eruption is due to a combination of mechanisms, including those 
mentioned above. 
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Fig. 4. The height-time profiles for the two eruptions of the sim- 
ulation. The curve for the first eruption is solid and the second is 
dashed. 



The height-time profiles of the erupting ropes are displayed 
in Figure [4] These are estimated by tracking the O-point of the 
magnetic field in the y = plane. The O-points also carry a 'cir- 
cular' profile of dense plasma so we are confident that these ac- 
curately represent the erupting flux ropes. This is demonstrated 
in Figure[5] which shows the second erupting flux rope, with cen- 
tre (x,z) ~ (0,55), at t = 132 by displaying the plasma density 
in the y — plane. Field line arrows are also shown to indicate 
the direction of the magnetic field in that plane. 

As mentioned earlier, the first flux rope forms at t « 80. It 
rises slowly until the eruption at t » 91. Here the gradient of the 
curve, in Figure |4] changes from (in non-dimensionalized units) 
0.3 to 2. At t = 115 the centre of the rope has reached z = 70 and 
the gradient of the curve has increased to 2.8. We stop tracking 
the rope at this time since this is just before it comes into contact 
with the top boundary of the computational box. At t « 125 the 
second rope, that has formed, is ejected upwards. Again, the rope 
is only tracked to just above z — 70 since beyond this the upper 
boundary begins to interfere with the rope's ascent. 

The simulation is ended after the second eruption due to 
computational constraints. The plasma from the first two erup- 
tions hits the top of the computational box, falls back down and 
interferes with the system. In theory, however, as long as the 
emergence process transports enough flux into the atmosphere 
and drives the necessary shearing, there should be more erup- 
tions like the ones described. 

The ropes are carried upwards by a Lorentz force. Looking 
at the vertical forces, (j x B) ; - dp/dz - pg > at the height of 
a flux tube during its rise. Figure [6] illustrates this by displaying 
the (non-dimensionalized) vertical forces as a function of height 
at (x,y) = (0,0) for the second flux rope eruption at t — 130. 
A positive upward force exists at the height of the flux rope and 
moves upwards in time with it. Of all the forces, it is the Lorentz 
force that dominates and is ultimately responsible for the rise of 
the ropes. 



4. Summary 

In this paper we have demonstrated that multiple CME-like 
eruptions are possible from a toroidal flux tube emerging into 
a magnetized corona. This com bin es and builds on the wor k 
of lMacTaggart & Hoodl (f2009bl) and lArchontis & TorEkl (12008). 
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Fig. 5. The second flux tube eruption of the simulation at t — 
132. The colour map shows logp in the y — plane with field 
line arrows to show the direction (but not the magnitude) of the 
magnetic field. Dense plasma can clearly be seen to be carried 
upwards into the corona with the erupting flux tube. 
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Fig. 6. The vertical forces carrying the flux rope upwards at t = 
130. The cut is taken at (x,y) = (0,0). Key: All forces (solid), 
(j x B) z (dash), -dp/dz (dot-dash). 



For the present study we consider a corona with a field that is 
almost antiparallel to the field of the emerging tube. External 
reconnection at the apex of the emerging arcade weakens the 
tension of the coronal field. With the expansion of the arcade, 
this reconnection becomes faster through time. Shearing, which 
occurs as part of the emergence process, induces internal recon- 
nection in the arcade and produces a flux rope. The continued 
emergence in combination with removal of the overlying coro- 
nal field eventually results in the expulsion of the flux rope. The 
mechanism for this eruption is similar to that of the breakout 
model. One important difference, however, is that the external re- 
connection in this model does not take place at a null point. After 
the first eruption, continued emergence and, therefore, shearing 
results in the formation of a second flux rope. This also erupts 
but the trigger mechanism cannot be directly linked to the break- 
out model, as with the first eruption. Due to the reconnection of 
the first eruption with the corona, a weakened coronal field ex- 
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ists above the second rope when it forms. Possible trigger mech- 
anisms, such as runaway reconnection and the torus instability, 
have been suggested. However, it is possible that the trigger for 
the second eruption is a combination of such mechanisms. 

Acknowledgements. DM and AWH would like to thank the referee, Tibor Torok, 
whose comments improved this research note. DM acknowledges financial as- 
sistance from STFC. The computational work for this paper was carried out on 
the joint STFC and SFC (SRIF) funded cluster at the University of St Andrews. 
DM and AWH acknowledge financial support from the European Commission 
through the SOLAIRE Network (MTRN-CT-2006-035484). 



References 

Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485 
Arber, T. D., Longbottom, A. W., Gerrard, C. L., & Milne, A. M. 2001, Journal 

of Computational Physics, 171, 151 
Archontis, V., Moreno-Insertis, F, Galsgaard, K., & Hood, A. W. 2005, ApJ, 

635, 1299 

Archontis, V. 2008, J. Geophys. Res., 1 13, Issue A3 
Archontis, V. & Torok, T. 2008, A&A, 492, L35 
Bateman, G. 1978, MHD Instabilities (MIT Press) 
DeVore, C. R. & Antiochos, S. K. 2008, ApJ, 680, 740 
Fan, Y. 2009, ApJ, 697, 1529 
Forbes, T. G. 2000, J. Geophys. Res., 105, 23153 

Galsgaard, K., Archontis, V., Moreno-Insertis, R, & Hood, A. W. 2007, ApJ, 
666,516 

Hood, A.W., Archontis, V., Galsgaard, K G. & Moreno-Insertis, F. 2009, A&A, 
503, 999 

Hornig, G. & Priest, E.R. 2003, Physics of Plasmas, 10, 2712 

Kliem, B. & Torok, T. 2006, Phys. Rev. Lett., 96, 255002 

Lynch, B. J., Antiochos, S. K., DeVore, C. R., Luhmann, J. G., & Zurbuchen, 

T. H. 2008, ApJ, 683, 1 192 
Mackay, D.H. & van Ballegooijen, AA. 2006, ApJ, 641, 577 
Maclean, R.C., Parnell, C.E., Galsgaard, K. 2009, Sol. Phys., accepted 
MacTaggart, D. & Hood, A.W. 2009a, A&A, 501, 761 
MacTaggart, D. & Hood, A.W. 2009b, A&A, accepted 

Manchester, W., IV, Gombosi, T., DeZeeuw, D. & Fan, Y. 2004, ApJ, 610, 588 
Schindler, K. 2008, Physics of Space Plasma Activity (Cambridge University 
Press) 

Soenen, A., Zucarello, F.P., Jacobs, C, Poedts, S., Keppens, R. & van der Hoist, 

B.2009, A&A, 501, 1123 
Torok, T. & Kliem, B. 2007, Astron. Nachr, 328, 743 
Vrsnak, B. 2008, Ann. Geophys., 26, 3089 



